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Abstract 

A distinct peak is observed in the Tr^ip' invariant mass distribution near 4.43 GeV in B -> Ktt^' 
decays. A fit using a Breit-Wigner resonance shape yields a peak mass and width of M = 4433 ± 
4 (stat) ± 2 (syst) MeV and T = 45^® ( stat )-i3 ( s y st ) MeV - The Product branching fraction is 
determined to be B(B° -► K^Z ± (U30)) x £(£±(4430) -► ttV) = (4.1 ± l.O(stat) ± 1.4(syst)) x 
10 -5 , where Z =t (4430) is used to denote the observed structure. The statistical significance of 
the observed peak is 6.5a. These results are obtained from a 605 fb -1 data sample that contains 
657 million BB pairs collected near the T(4S I ) resonance with the Belle detector at the KEKB 
asymmetric energy e + e~ collider. 

PACS numbers: 14.40. Gx, 12.39.Mk, 13.25.Hw 
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An unexpected bonus from the 5-factory experiments has been the discovery of a large 
number of charmonium-like meson states p], 0, 0, 0, 0, 0, 0, 0, 0| . While two of these have 
been identified as likely candidates for "missing" charmonium states, namely the rj c [l| and 
the x'c2 0; others have properties that are at odds with expectations of the charmonium 
model. These latter include the X(3872) Q and V(4260) fl, which are seen to decay to 
7i + 7i-J/ip; the X(3940) @, seen in D*D; the F(3940), seen in luJ/^j Q; and the V(4325) Q, 
seen in 7r + 7r~il)'. Recently Belle reported a second 7i + 7i^J/ip mass enhancement below the 
V(4260) [8] and has shown that the V(4325) — > n^n^tp' signal, reported by BaBar to have 
a width of T = 172 ± 33 MeV, is better fitted with two narrower peaks, one at 4361 MeV 
with r = 70 ± 20 MeV and a second at 4664 MeV, with T = 40 ± 17 MeV Q. 

Proposed assignments for these states have included: multiquark states, either of the 
(cq, cq) "molecular" type [lOj or [cq, cq] diquark-antidiquark type [rjj (here c represents a 
charmed quark and q either a u-, d- or s-quark); hybrid cc-gluon mesons [12j; or other missing 
charmonium states where the masses predicted by potential models are drastically modified 
by nearby D^D^ thresholds 13j, [l4j. A characteristic that clearly distinguishes multiquark 
states from hybrids or charmonia is the possibility to have charmonium-like mesons with 
non-zero charge (e.g. [cued]), strangeness ([edes]) or both ([cues]) [151 ] . These considerations 
motivated a search for charmonium-like mesons with non-zero electric charge. 

Here we report the observation of a relatively narrow peak in the ir + %p' invariant mass 
distribution produced in exclusive B — > Kn + ip' decays [16]. The results are based on an 
analysis of a 657 million BB event sample collected in the Belle detector operating at the 
KEKB asymmetric energy e + e~ collider. The data were accumulated at a center-of-mass 
system (cms) energy of yfs = 10.58 GeV, corresponding to the mass of the T(45") resonance. 
KEKB is described in detail in Ref. 17| . 



The Belle detector, described in Ref. [18], is a large-solid-angle magnetic spectrometer 
that consists of a silicon vertex detector, a 50-layer cylindrical drift chamber, an array of 
aerogel threshold Cherenkov counters, a barrel-like arrangement of time-of-flight scintillation 
counters, and an electromagnetic calorimeter comprised of CsI(Tl) crystals located inside 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K\ mesons and to identify muons. 

We select events of the type B — ► Ktt + iP', where the ip' decays either to £ + £~ or 7r + 7r~ J/ip 
with J/xjj — > £ + £~ {£ = e or fi). Both charged and neutral (K s — > 7r + 7r~) kaons are used. 
Charged tracks other than K s — ► 7r + 7r~ secondaries are required to originate from the 
beam-beam interaction point. The charged kaon, lepton and pion selection requirements are 
described in Ref. 0; those for neutral kaons are described in Ref. 0. 

For ip'(J/ip) — > £ + £~ candidates we require the invariant mass of the lepton pair to be 
within 20 MeV of the ijj'(J/ijj) mass. For ip'(J/ip) — > e + e~ candidates, we include photons 
that are within 50 mrad of the e + or e~ tracks in the invariant mass calculation. For 
if)' — > 7t + 7t~J/iI) candidates, we require the 7r + 7r~ invariant mass to be greater than 0.44 GeV 
and \M(tc + ti-£ + £-) - M(£ + £~) - 0.589 GeV| < 0.0076 GeV, which is ±2.5(7, where a is the 
rms resolution. 

We suppress continuum e + e~ — > qq events, where q = u, d, s or c, by requiring R 2 < 0.4, 
where R 2 is the second normalized Fox- Wolfram event-shape moment [l9[ . We also req uire 
| cos#b| < 0.9, where 9b is the angle between the B meson and e + beam directions [20]. 



We identify B mesons using the beam-constrained mass M bc = J E^ CSLm — p 2 B and the 
energy difference AE = E heam — E B , where -Ebeam is the cms beam energy, p B is the vector 
sum of the cms momenta of the B meson decay products and Eb is their cms energy sum. 
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We select events with \M^ C — m^l < 0.0071 GeV (tub = 5.279 GeV, is the world average 
5-meson mass 2JJ) and \AE\ < 0.034 GeV, which are ±2.5cr windows around the nominal 
peak values. 

The invariant mass of the selected B — > Knif)' candidate tracks is kinematically con- 
strained to equal rag. This improves the if)' — > £ + £~ (J/if> — > £ + £~) mass resolution to 
a = 4.4 MeV (5.3 MeV). We require M(£ + £~) computed with the fitted lepton four- vectors 
to be within ±2.5cr of m^> (m^), the world average if>' (J/if>) mass (2lj . 

For the ip' -> £ + £~ mode we compute M(mf)') as M(n£ + £-) - M(£ + £~) + m f ; for if)' -> 
7r + 7r~ J /if; decays, we use M(7rif>') = M(it'K + 'K~J/if))—M('K + 'K~J/if)) J r'm^i. Simulations of the 
two if)' decay modes indicate that the experimental resolution for M(n + if)') is a ~ 2.5 MeV 
for both modes. 




FIG. 1: The M 2 (Ktt) (horizontal) vs. M 2 (vn//) (vertical) Dalitz-plot distribution for B° -> 
K~Tr + ip' candidate events. 

Figure [1] shows a Dalitz plot of M 2 (Ktt + ) (horizontal) vs. M 2 (iT + if)') (vertical) for the 
B — > KiT + if)' candidate events. Here, a distinct band at M^ n ~ 0.8 GeV 2 , corresponding 
to B -> K*(890)ip'; K*(890) -> Kn, is evident. In addition, there are signs of a ^(1430) 
signal near M^ n = 2.0 GeV 2 . The B — > IT* (890) V>' events are used to calibrate the M bc and 
A£7 peak positions and widths. 

Some clustering of events in a horizontal band is evident in the upper half of the Dalitz plot 
near M 2 (7n/> / ) ~ 20 GeV 2 . To study these events with the effects of the known Kir resonant 
states minimized, we restrict our analysis to the events with \M(Ktt) —ttlk* (89o) \ > 0.1 GeV 
and \M(Kir) — mx*(i43o)| — 0.1 GeV. In the following, we refer to this requirement as the 
K* veto. 

The open histogram in Fig. [2] shows the M(n + if)') distribution for selected events with 
the K* veto applied. The bin width is 10 MeV. The shaded histogram shows the scaled 
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distribution from AE sidebands (\AE ± 0.070| < 0.034 GeV). Here a strong enhancement 
is evident near M(irip') ~4.43 GeV. 
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FIG. 2: The M(tt + ^') distribution for events in the M^ C -AE signal region and with the K* veto 
applied. The shaded histogram show the scaled results from the AE sideband. The solid curves 
show the results of the fit described in the text. 

We perform a binned maximum-likelihood fit to the M(jnp') invariant mass distribution 
using a relativistic S'-wave Breit-Wigner (BW) function to model the peak plus a smooth 
phase-space-like function / CO nt(M), where / con t(M) = A/" con t<f (Q 1/2 + ^4iQ 3/2 + A 2 Q 5/2 ). 
Here q* is the momentum of the ir + in the nip' rest frame and Q = M max — M, where 
M max = 4.78 GeV is the maximum M(jnp') value possible for B — > Knip' decay. The 
normalization 7V con t and two shape parameters A\ and A2 are free parameters in the fit. 
This form for f CO nt{M) is chosen because it mimics two-body phase-space behavior at the 
lower and upper mass boundaries. (Since the M(jnp') distribution for the non-peaking B- 
decay events and the AE sideband events have a similar shape, we represent them both 
with a single function.) 

The results of the fit, shown as smooth curves in Fig. [2], are tabulated in Table [B The 
fit quality is x 2 = 80-2 for 94 degrees of freedom. The significance of the peak, determined 
from the change in log likelihood when the signal and its associated degrees of freedom are 
removed from the fit, is 6.5cr. 



TABLE I: Results of the fit shown in Fig. H 





-V C ont 


BW Mass (GeV) 


T (GeV) 


121 ± 30 


766 ± 39 


4.433 ± 0.004 


u.u4o_ 013 
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We fit M(nip') distributions for various subsets of the data. The results are summarized 
in Table HH 

There are significant (i.e. significance of 4.5<r or more) signals with consistent mass values 
in both the if}' — > 7i + 7i^J/ip and if}' — > £ + £~ subsamples. However, the width of the peak 
in the if>' — > £ + £~ subsample is substantially wider than that for the n + ir~J/ifj subsample. 
Fitting the two measured widths to a common value gives a x 2 — 4.8 for 1 degree of freedom. 
The corresponding confidence level is ~ 3%. 

The fitted values for the signal yields are highly correlated with the widths. To compare 
the yields in each subchannel, we refit the distributions using a width that is fixed at the 
T = 0.045 GeV value determined from the common fit. The ratio of "constrained" signal 
yields for the 7i + ii~J/ip and £ + £~ subsamples (see Table ITT]) is 1.09±0.35, in good agreement 
with the MC-determined acceptance ratio of 1.23. 



TABLE II: Results of fits to different subsamples of the data. 



Subset 


Mass 


Width 


Signif. 


Constr. yield 




(GeV) 


(GeV) 


(*) 


(r = 0.045GeV) 


Tr + 7T~J/lp 


4.435 ± 0.004 


0-026^ 


4.5 


64 ±15 


e+r 


4.435 ±0.010 




4.7 


59 ±13 




4.430 ± 0.009 


0.056±g : gg 


3.5 


41 ±12 


n + n~ 


4.434 ± 0.004 


o.o38i8:85i 


5.2 


80 ±16 


K^ifj' 


4.434 ± 0.005 


o.o48±8:ga 


6.0 


102 ± 18 




4.430 ± 0.009 


0.048-fixed 


2.0 


19 ±8 


K* veto 


4.437 ± 0.005 


0.0631833 


7.1 


170 ± 26 



Table [IT] also shows the results from dividing the data sample into £ + £~ = e + e~ and 
yU + /i _ , and charged kaon and Kg — > 7r + 7r~ subsets. We see signals in both the e + e~ (3.5c) 
and /^ + yU~ modes (5.2a) with consistent mass and width values, and with constrained yields 
that are consistent with the expected e + e~/fi + fi~ acceptance ratio of 0.61. There are too 
few events in the K s sample to enable a stable fit with yield, mass and width all allowed to 
vary. With the width fixed at the value found for the charged kaon sample, the fit returns a 
19 ± 8 event signal with 2.0<r significance and a consistent mass value. The observed signal 
yield in the K s sample agrees with expectations 22| based on scaling the charged kaon 
signal by the Kg/K^ acceptance ratio (0.19). 

The last row of Table HT1 shows the results of a fit to the M(7n//) distribution for the case 
where the K* veto is replaced by a less stringent requirement that only eliminates the core 
of the K*(890) peak: \M(Kn) — mx*(890)| > 0.05 GeV. Here the observed signal increases 
and its statistical significance improves to 7.1er. 

The M(Kir) distribution for events within ±0.03 GeV of the peak at 4.43 GeV is shown 
in Fig. [3J Here the K* veto, which excludes the regions indicated by the double-sided arrows 
in the figure, has been removed. The shaded histogram is the scaled AE sideband data. 
Aside from the if* (890) resonance events, which are removed by the K* veto, no dramatic 
features are evident. 

We considered the possibility that interference between S-, P- and D-waves in the Ktt 
system might produce a structure similar to that which is observed. (There are F-wave 
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FIG. 3: The M(Ktt) distribution for events within ±0.03 GeV of the 4.43 GeV peak. Here the K* 
veto, which excludes the regions indicated by double-sided arrows, has been removed. The shaded 
histogram indicates the AE'-sideband-determined non-1? background. 

and higher Kit resonances listed in the PDG tables, however, even the lowest mass F-wave 
entry, the ^3 (1780), is not kinematically accessible in B — > Knip' decay.) We find that 
with only these three partial waves, it is not possible to produce a nip' invariant mass peak 
near 4.43 GeV that is as narrow as the one we see without other, even more dramatic, 
accompanying structures. 

We applied the same analysis to large MC samples of generic B meson decays and found 
no evidence of peaking in the 7iip' invariant mass distribution. 

The product branching fraction is determined using MC-computed acceptance values and 
world average values for ip' and J/ip branching fractions (2l| . For this calculation, we only 
use the signal yield from the B° — * K T TT ± ifj' decay sample. The resulting product branching 
fraction is 

B{B° -> F~Z + (4430)) x £(Z+(4430) -> tt + <//) 

= (4.1 ±1.0 ±1.4) x 10 -5 , (1) 

where Z + (4430) is used to denote the observed structure, the first error is statistical and 
the second error is systematic (discussed below). 

The values of the ip' mass determined using events in the data agree within ±1 MeV with 
the world average value for both the ip' — ► 7t + tt~ J/ip and ip' — > £ + £~ decay modes. We find 
less than 1 MeV variation in the peak mass value for different fitting functions. Fits that 
include possible interference between the BW signal and the non-resonant irip' continuum 
produce at most a 1.2 MeV shift in the fitted mass value. We assign a ±2 MeV systematic 
error to the mass determination. 

The systematic uncertainty on the width is mostly due to the uncertain effects of back- 
ground fluctuations feeding into the fitted signal. We estimate the level of this effect from 
the range in width values determined from different subsets of the data to be +29 MeV 
and —10 MeV. Changes in the parameterization of f CO nt{M) and variations in the range of 
M(irip') values included in the fit produce ±8 MeV changes in the width; using different 
BW forms produce ±4 MeV width changes. Adding these sources in quadrature results in 
a total systematic error on the width of ±f z MeV. 

The largest systematic error on the product branching fraction measurement is due to 
the correlation between the fitted signal yield and the peak width. A +30 MeV (—13 MeV) 
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change in the width produces a +30% (—17%) change in signal yield. Changes in the 
parameterization of f CO nt(M) produce +11% and —18% variations in signal yield. Other 
systematic errors are smaller. These include: possible interference with the nip' continuum; 
the choice of BW signal function; uncertainties in the acceptance calculation; uncertainties 
in the tracking and particle identification efficiencies; errors on the world-average if)' decay 
branching fractions; MC statistics; and the error on the number of BB mesons in the sample. 
Combining these errors in quadrature gives a systematic error on the product branching 
fraction of 35%. 

In summary, a study of B — > Kir + ifj' decays reveals a peak in the 7r +, 0' invariant 
mass spectrum at M — (4433 ± 4(stat) ± 2(syst)) MeV. The measured width, T = 
(45l]J(stat)^3(syst)) MeV, is too narrow to be caused by interference effects in the Ktc 
channel. The statistical significance of the observed peak is 6.5<r. 

There have been a number of anomalous charmonium-like meson candidates reported in 
the literature 0, 0, 1, I, H, @, B @, B- The structure reported here is unique in that it is 
the first candidate to have a non-zero electric charge. 
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MOEHRD, KOSEF and KRF (Korea); KBN (Poland); MES and RFAAE (Russia); ARRS 
(Slovenia); SNSF (Switzerland); NSC and MOE (Taiwan); and DOE (USA). 



[1] S.-K. Choi et al. (Belle Collaboration), Phys. Rev. Lett. 89, 102001 (2002). 
[2] S. Uehara et al. (Belle Collaboration), Phys. Rev. Lett. 96, 082003 (2005). 
[3] S.-K. Choi et al. (Belle Collaboration), Phys. Rev. Lett. 91, 262001 (2003). 
[4] B. Aubert et al. (BaBar Collaboration), Phys. Rev. Lett. 95, 142001 (2005). 
[5] K. Abe et al. (Belle Collaboration), Phys. Rev. Lett. 98, 082001 (2007). 
[6] S.-K. Choi et al. (Belle Collaboration), Phys. Rev. Lett. 94, 182002 (2005). 
[7] B. Aubert et al. (BaBar Collaboration), hep-ex/0610057, submitted to Physical Review Let- 
ters. 

[8] C.Z. Yuan et al. (Belle Collaboration), arXiv:0707.2541, submitted to Physical Review Letters. 
[9] X.L. Wang et al. (Belle Collaboration), Phys. Rev. Lett. 99, 142002 (2007). 
[10] See, for example, M.B. Voloshin and L.B. Okun, JETP Lett. 23, 333 (1976); M. Bander, 

G.L. Shaw and P. Thomas, Phys. Rev. Lett. 36, 695 (1977); A. De Rujula, H. Georgi and 

S.L. Glashow, Phys. Rev. Lett. 38, 317 (1977); N.A. Tornqvist, Z. Phys. C 61, 525 (1994); 

A.V. Manohar and M.B. Wise, Nucl. Phys. B 339, 17 (1993); N.A. Tornqvist, Phys. Lett. 

B590, 209 (2004); F.E. Close and PR. Page, Phys. Lett. B 578, 119 (2003); C.-Y. Wong, 

Phys. Rev. C 69, 055202 (2004); E. Braaten and M. Kusunoki, Phys. Rev. D 69, 114012 

2004); E.S. Swanson, Phys. Lett. B588, 189 (2004). 
[11] T.-W. Chiu and T.-H. Hsieh, Phys. Lett. B646, 95 (2007). I. Bigi et al, Phys. Rev. D 72, 

114016 (2005); L. Maiani et al, Phys. Rev. D 72, 031502 (2005). 
[12] D. Horn and J. Mandula, Phys. Rev. D 17. 898 (1978); F.E. Close, Phys. Lett. B 342, 369 

(1995); C. McNeile, C. Michael and P. Pennanen (UKQCD Collaboration), Phys. Rev. D 65, 

094505 (2002). 



9 



[13] T. Barnes and S. Godfrey, Phys. Rev. D 69, 054008 (2004); E.J. Eichten, K. Lane and 

C. Quigg, Phys. Rev. D 69, 094019 (2004); M.B. Voloshin, Phys. Lett. B579, 316 (2004). 
[14] C. Meng and K.-T. Chao, Phys. Rev. D 75, 114002 (2007). 
[15] T.-W. Chiu and T.H. Hsieh, Phys. Rev. D 73, 111503 (2006). 

[16] In this letter the inclusion of the charge conjugate mode is always implied. The tp' is sometimes 
denoted as ip(2S). 

[17] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 479, 117 (2002), and other papers 

included in this volume. 
[18] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A479, 117 (2002). 
[19] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[20] The 9b distribution has a sin 2 #b dependence for e + e~ — > BB and is flat for continuum 

background events. 
[21] W.M. Yao et al. (Particle Data Group), J. Phys. G 33, 1 (2006). 
[22] This assumes equal branching fractions for B° and B + mesons. 



10 



